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The review of the nucleon spin structure functions problems is pre- 
sented. The perturbative QCD predictions for the small x behaviour of the 
nucleon spin structure functions is discussed. The role of the resummation 
of the In 2 \ jx terms is emphasised. Predictions for the nonsinglet structure 
function g\ in case of a flat as well as a dynamical input are given. The 
so called 'spin crisis' in the context of both theoretical and experimental 
future hopes are also briefly discussed. 
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1. Introduction 

Since 1988, when the famous EMC experiment [1] provided surprising 
results, the polarised deep inelastic lepton-nucleon scattering (DIS) became 
very interesting from experimental as well as theoretical point of view. This 
experiment, in which longitudinally polarised muons scattered on longitu- 
dinally polarised protons, brought the conclusion that quarks are carrying 
only a small part of the proton spin projection in the polarised proton. 
This result called as 'spin crisis' is still a challenge for theoretical and ex- 
perimental research. The main questions should be answered are: how the 
nucleon spin is distributed among its constituents: quarks and gluons and 
how the dynamics of these constituent interactions depend on spin. Solu- 
tions of these problems may be found within perturbative QCD because they 
involve hard and semihard (short-distance) processes. From experimental 
point of view there are many projects, which provide or will provide a mech- 
anism to probe the spin properties of nucleon and should be a crucial test 
of QCD. The most important of these projects are experiments in SLAC, 
DESY (HERMES, HERA), CERN (SMC) (with polarised proton, deuteron 
and He 3 targets), CERN (COMPASS) (with polarised muon beams and 
longitudinally polarised hydrogen (NH 3 ) and deuteron (Li 6 D) targets. The 
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main goal of all these experiments is to measure the nucleon spin structure 
functions g\(x,Q 2 ) and g2(x,Q 2 ). It can be done by measurement of the 
cross section asymmetry, where two considered cases correspond to antipar- 
allel or parallel spin orientation of the longitudinally polarised lepton (//, e) 
and nucleon (p, n, d) respectively. Recently the experiment data has allowed 
to investigate the nucleon spin structure in the large range of the kinemat- 
ical variables: Bjorken x and Q 2 . The most interesting, both theoretically 
and phenomenologically, is the region of small x. Theoretical understanding 
of the small x (x ~ 10~ 3 and less) behaviour of the polarised nucleon struc- 
ture function enables the correct estimation of Ti momenta in sum rules. It 
is very important because present experimental data do not cover however 
the whole very small x region and the only way (at present) to know the 
nucleon spin structure completely is extrapolation of large and medium x 
results into the small x region through the theoretical QCD analysis. From 
the other side, future polarised experiments in HERA [2] will enable spin 
DIS investigations in the very small x region: x ~ 1CT 4 and less. Then 
theoretical predictions would be verified by the experiment. These future 
spin experiments would be a crucial test of theoretical analysis. Descrip- 
tion of the nucleon spin structure function g\ within perturbative QCD for 
small x can be done in different frames (in LO, NLO, In 1/x, In 2 1/x etc. 
approximations) giving different results for g± in this region. Thus the fu- 
ture comparison of theoretical and experimental results could be definitive. 
In the next section we shall briefly remind the sum rules and the 'spin cri- 
sis' problem. In section 3 we shall discuss the polarised structure functions 
of nucleon in the small Bjorken x region. We shall emphasise the In 2 1 jx 
resummation which is significant in this region. In point 4 the nonsinglet 
g^ s (x, Q 2 ) predictions are presented. We show LO and unified LO+ln 2 1/x 
resummation results in case of a flat (nondynamical) and a dynamical input 
parametrisation as well. We compare our results with recent SMC data. Fi- 
nally in conclusions we shall briefly discuss future experimental hopes and 
possible scenario of solving the spin crisis problem. 



There are four basic nucleon structure functions: F±, F2 for a spin inde- 
pendent case and g±, 52 for a polarised one, which characterise the pointlike 
interaction between virtual hard photon (Compton scattering) of Q 2 A 2 
and hadron constituents - partons in the deep inelastic lepton - hadron 
scattering. In the parton model Fi , F2 and g\ have a very simple form and 
interpretation. Thus 
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where is a charge of the i-flavour quark, ft+(a?) (ft_(x)) is the density 
distribution function of the i-quark with the spin parallel (antiparallel) to 
the parent nucleon. In the Bjorken limit 

F 2 (x) = 2xF 1 (x) (2.2) 



9i{x) = \ E e ? A <?^) ( 2 -3) 

i=u,d,s,.. 

Aft(x) =q i+ (x) -qi-(x) (2.4) 

Function gi(x,Q 2 ) is connected with the helicity of the nucleon (i.e. spin 
projection on the momentum direction). Thus the integral 

l 

<A ft ) = J A qi (x)dx (2.5) 
o 

is simply a part of the nucleon helicity, carried by a quark of i-flavour 
(i=u,d,s,..). The second spin dependent structure function g2{x, Q 2 ), related 
to the transverse spin polarisation of the nucleon has no simple meaning in 
the parton model. The main goal in the deep inelastic lepton - nucleon 
scattering experiments with polarised both the lepton and the nucleon par- 
ticles is to find the spin dependent structure function of nucleon gi(x,Q 2 ). 
This measurement of g±(x, Q 2 ) provides the knowledge how the spin of the 
nucleon is distributed among the partons: valence quarks u v ,d v , sea quarks 
q sea and gluons g. The experimental measurement of g\(x : Q 2 ) is based on 
the measurement of the cross section asymmetry factor A [13]: 

A = ^_^±d (2 . 6) 

where a(++), a(-\ — ) correspond to the cases when the spins of longitu- 
dinally polarised lepton (jj, or e) and nucleon (p or n or d) are parallel or 
antiparallel. Finally, g\(x,Q 2 ) can be determined through the relation 

where 

Al = 01/2 - ° 3/2 (2.8) 

0"l/2 + &3/2 

and 

R = — (2.9) 
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The cross section a 3 / 2 an d &1/2 correspond to the helicities 3/2 and 1/2 of 
the absorbed virtual photon. &l and ut are the absorption cross sections 
of longitudinal and transverse virtual photons in the polarised scattering. 
Through polarised DIS experiments with different targets: proton, neutron 
or deuteron ones, it is possible to combine the g\ results for these different 
cases and hence to find out the spin dependent distribution functions of par- 
tons (Au v , Ad v , Aq sea , Ag) in the nucleon. The measurement of g\ function 
enables also verification of the sum rules, which play a very important role 
as a test of QCD. The most significant of them are Bjorken and Ellis-Jaffe 
sum rules [7]. Both sum rules are related with the nucleon spin structure 
functions and their estimation requires the knowledge of first moments of 
gi for proton, neutron and deuteron T'{(Q 2 ), r?(Q 2 ), Tf(Q 2 ): 

1 

r\(Q 2 ) = J g\{x,Q 2 )dx (2.10) 


To test experimentally the Bjorken and Ellis-Jaffe sum rules it is neces- 
sary to know function g\(x,Q 2 ) for certain Q 2 and in the entire range of 
x: x 6 (0; 1). This data collection of g 1 for arbitrary values of Q 2 and x 
variables is however impossible because of technical constraint in the ex- 
periments. The broad ranges of x and Q 2 (where Q 2 = —q 2 , q is the four 
momentum transfer between lepton and hadron and W 2 = Q 2 (l/x — 1) is 
total CM energy squared) cannot be reached independently. The accessible 
at present kinematical region in the polarised fixed target HERMES ex- 
periment in HERA is 0.004 < x < 1, 0.2 < Q 2 < 20 GeV 2 [26]. It must 
be however emphasised, that HERA will extend the present region by two 
orders of magnitude both in x and Q 2 : x down to 6 • 10~ 5 and Q 2 up to 
2 • 10 4 GeV 2 . It will take place in the future experiment with polarised 
electron and proton beams [8]. Hence physicists hope for a new very inter- 
esting field of experimental investigations. Broader range of x and Q 2 in the 
future polarised DIS experiments in HERA will enable more precise knowl- 
edge of the proton spin structure and the test of QCD predictions for the 
polarised structure functions of the nucleon. The accessible region of x and 
Q 2 in future polarised experiments in HERA is roughly presented in Fig.l 
[8]. Perturbative QCD provides methods for theoretical analysis of struc- 
ture functions for Q 2 » A 2 (A ~ 200 MeV). When the coupling constant of 
strong interactions a s (Q 2 ) becomes large (a s ~ 1) perturbative calculations 
are not applicable. QCD as a theory with asymptotic freedom improves the 
naive parton model by permission for a weak Q 2 dependence of structure 
functions F±, F2, gi, 52- This scaling violation is caused by strong inter- 
actions between partons (quarks and gluons). The larger Q 2 (i.e. smaller 
distance between partons) the weaker interaction between partons. Pertur- 
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bative QCD enables Q 2 evolution of structure functions, where the values 
of the structure functions at some fixed scale Q 2 = Qq are taken from the 
experiment and used as an input in evolution equations. As it was already 
mentioned above, the most important test of QCD in the spin dependent 
DIS phenomena is verification of Bjorken (BSR) and Ellis-Jaffe (EJSR) sum 
rules. BSR is related to the isovector axial current SU(3) flavour symme- 
try in the nucleon f3 decay while EJSR results from the octet axial current 
SU(3) symmetry in the (3 decay of hyperons. In the parton model BSR and 
EJSR have forms: 

BSR : a 3 = g A (2.11) 

where 

l 

a 3 = J (An + Au - Ad - Ad)dx (2.12) 
o 

and gA ~ 1.257 is (3 decay structure constant for neutron; 



EJSR : a$ = 3F — D 



(2.13) 



G 



where 

1 1 

a 8 = J g° l ctet (x)dx = j{Au + An + Ad + Ad- 2As - 2As)dx (2.14) 
o o 

Aq(q = u,d,s,..) are polarised distribution functions of quarks and an- 
tiquaries (2.4) and F, D are axial coupling constants, estimated from the 
weak (3 decays of hyperons: 3F — Dm 0.579. Taking into account the QCD 
corrections one can obtain: 

ao,3,8 -»■ oo,3,8(l - corr(a)) (2.15) 

where corr(a) is just a QCD correction, calculated in a perturbative way 
[13]: 

corr(a) w (a s /it) + 3.58(a s /vr) 2 + 20.22(a s /vr) 3 + 130(a s /vr) 4 (2.16) 

a s (Q 2 ) is the running coupling constant of strong interaction and in LO 
approximation has a form 

a s (Q2 } = 12t[ ( } 

(33-2^)^^ 

where A = 232 MeV is a scale parameter of QCD. ao in (2.15) denotes the 
singlet part of axial current: 

l l 
a = J g s i n9le \x)dx = J AZdx (2.18) 


For three flavours of quarks (Nf = 3) AS is given by 

Aq s = AS = An + Au + Ad+ Ad+ As + As (2.19) 

Thus, taking into account QCD corrections (2.15) one can rewrite the BSR 
(2.11) in the form 

l 

r? - r? = f (gl - 9?)dx = l -g A {l - corr(a)) (2.20) 



Bjorken and Ellis- Jaffe sum rules imply the following expressions for first 
moments of spin structure functions g\ of the nucleon: 

l 

r^J g^dx = (g + | + |)(1 " corr(a)) - ^a s (Q*)(Ag{Q*)) 
o 

(2.21) 
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l 

T^J 9?(x)dx = (-g + H + |)(1 - corr(a)) - -^a s (Q 2 )(A 5 (Q 2 )> 
o 

(2.22) 

l 

It = (1 - \u D ) J gf(x)dx = [(g + - corr{a)) 

o 

-^a s (Q 2 )(A 5 (Q 2 ))](l-^) (2.23) 

where Ag is the spin dependent distribution function of gluons and (Ag) 
is defined similarly as for quarks (2.5). Parameter lod ~ 0.058 denotes 
probability that deuteron is in D state. QCD corrections a™ in formulae 
(2.20)-(2.23) enable us to compare theoretical perturbative QCD predictions 
for the structure function g\ and its first moment Y\ with the experimental 
results (2.7). Thus one can easily read from (2.20), (2.16) and (2.17) that 
BSR at Q 2 = 10 GeV 2 for A = 200 MeV up to a* correction gives 

THEORY : — T™ = 0.185 (2.24) 

This result is in a good agreement with the most recent SMC and SLAC 
measurements [13], [14], which yield for Q 2 = 10 GeV 2 the quantity 

EXPERIMENT: — T" = 0.195 ± 0.029 (2.25) 

It means that the Bjorken sum rule BSR is fulfilled within the experimental 
errors. For the Ellis-Jaffe sum rule EJSR however, there is a significant 
disagreement between theory and experiment. Theoretical value for EJSR 
from the simplest version of (2.21), where the gluon contribution is neglected 
and the strange sea quarks are unpolarised (As = 0) gives [13] 

THEORY : T\{Q 2 = 10.7 GeV 2 ) = 0.171 (2.26) 

while the EMC result [1] is 

EXPERIMENT : T\{Q 2 = 10.7 GeV 2 ) = 

= 0.126 ± 0.010{statist) ± 0.015(system) (2.27) 

Comparing (2.26) and (2.27) one can see, that the EJSR prediction dis- 
agrees with experimental data by 2.6 standard deviations. From the latest 
experimental analyses [11], [12], [14] it can be read that BSR is always ful- 
filled within the experimental errors while EJSR is broken at the level of 
about 3 standard deviations. This fact is confirmed by many polarised DIS 
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experiments: SLAC in CERN with the polarised electron beam scattered 
on polarised proton (p), neutron (n) or deuteron (d) targets, SMC at CERN 
with muons fi + and p,d targets and HERMES at DESY experiments with 
positons e + and n,p,d targets. The characteristic of kinematic variables for 
experimental groups EMC, SLAC, SMC and HERMES is presented below 
in Tab.l [12]. The experimental data of g\, Ti and (Atfe) (2.5) are widely 



TABLE I 



Experiment 


Beam 


Lepton 
energy 
E (GeV) 


Smallest x 
Q 2 >1 
(GeV 2 ) 


Average 

Q 2 

(GeV 2 ) 


Nucleon 
target 


EMC 




100-200 


0.01 


10.7 


P 


SMC 




100-190 


0.003 


10 


p,d 


SLAC E-142 


e 


19-25 


0.03 


2 


n 


SLAC E-143 


e~ 


10-29 


0.029 


3 


p,d 


SLAC E-154 


e 


48.3 


0.014 


5 


n 


SLAC E-155 


e~ 


48.3 


0.014 


5 


p,d 


HERMES 


e+ 


27.5 


0.023 


2.3 


n,p,d 



Table 1. The kinematic variables in DIS for different experimental groups 



reviewed e.g. in [6], [16]. The main problem nowadays is to find out how the 
spin of the proton is distributed among partons: valence quarks, sea quarks 
and gluons. The part of the proton spin carried by the parton p, where p 
denotes u vah d vat , u sea , d sea , s sea , antiquarks, g is given by (2.28). 

l 

(Ap) = J Apdx (2.28) 
o 

and 



hence 



Aq = Aq val + Aq sea 
Aq sea = Aq sea = Aq 

Aq val = Aq - Aq 



(2.29) 
(2.30) 

(2.31) 



The method of the extraction Ap and then the important quantity (Ap) 
from experimental data is as follows: 

1. Polarised distribution functions Ap(x,Q 2 ) are parametrised at the given 
low scale Qq (e.g. Qq = 1 GeV 2 ) in the form: 



Ap(x,Q 2 ) = Nr] f x a f (1 - xff (I + -f f x s f) (2.32) 
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where 

1 

r, f = J Ap(x,Q 2 )dx (2.33) 
o 

N is a normalisation factor and a/, (3f, 7/, 5f are parameters. 

2. From QCD evolution equations (see Appendix A) (A.1)-(A.3) one can 
get in NLO approximation Ap(x,Q 2 ) functions for (x,Q 2 ) set, available in 
experiments. 

3. Using the spin dependent distribution functions Ap, one can calculate 
the spin structure function g\ (x, Q 2 ) from the formula: 

1 

gi (x,Q 2 )= l -{e 2 ) J ^-[Cl(x/y,a s (Q 2 ))A^y,Q 2 ) + 

X 

2N f C^x/y,a s (Q 2 ))Ag(y,Q 2 ) + C q NS (x/y,a s (Q 2 ))Aq N s(.y,Q 2 )} (2.34) 
where 

! N f 

(« 2 ) = ^E4 (2-35) 

= AS (2.19) is singlet and Aqjys is the nonsinglet polarised distribution 
function of quarks: 

N f 2 

Aq NS = E(yjy " !)( A % + (2-36) 

j=i > ' 

4. Now one can compare the calculated g\ (2.34) with experimental data for 
gi (2.7) and fit free parameters (a/, (3j, 7/, 5f) of input parametrisations 
(2.32) in such a way to minimise the x f° r an experimental points (x, Q 2 ). 
In this way it was found, as it was already mentioned above, that all ex- 
periments confirm the validity of the Bjorken sum rule while the Ellis-Jaffe 
sum rule is violated. After extraction of spin dependent distribution func- 
tions Ap from experimental data it was possible to determine the value 
J AE(s, Q 2 ))dx, which is a part of the proton spin carried by quarks. The 
result was surprising because it turned out that quarks carry only a very 
small part of the proton spin: 

1 

J AE(ar, Q 2 )dx w 0.17 ± 0.17 (2.37) 
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and the strange quarks carry much larger part of the proton spin that it 
follows from the simple parton model: 

l 

J As(x,Q 2 )dx « -0.14 (2.38) 
o 

This result known since EMC experiments [1] as a "spin crisis" problem 
should be correctly solved within QCD. Theoretical analysis should provide 
a proper interpretation of the EMC results. The spin problem is in fact 
not the "spin crisis" but the problem of understanding how the nucleon 
spin is composed of parton spins? Disagreement of experimental data with 
theoretical predictions, emerging in EJSR violation, shows that the simple 
quark model is definitely inadequate to the proper description of the nucleon 
spin structure. The "spin crisis" has arisen because experimental results did 
not confirm the naive quark model in which the proton spin is carried only 
by valence quarks i.e.: 

l l 

J AZdx = J (Au val + Ad val )dx = 1 (2.39) 



In this model it was assumed that strange sea quarks as well as gluons do 
not or almost do not contribute to the nucleon spin. Using this assumption, 
the input distribution functions (2.32) were parametrised and then used 
in the fitting of experimental data. So to remove the spin problem it is 
necessary to revise the previous knowledge about the spin distribution in 
the nucleon. It would be helpful with this aim to consider the following 
facts and assumptions [6]: 

1. The anomalous dimension of the first moment of gluons is equal to -1: 
l 

J Ag(x,Q 2 )dx ~ InQ 2 for Q 2 -> oo (2.40) 
o 

It means that 

l 

lim [a s (Q 2 ) [ Ag(x,Q 2 )dx] = const (2.41) 



and gluons give a nonvanishing correction of order a s to EJSR (2.21)-(2.23), 
when Q 2 — > oo (opposed to the first moment of quarks, which the anomalous 
dimension is zero) . This gluon anomalous dimension via evolution equations 
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(A.1)-(A.3) changes also the spin dependent distribution functions of quarks. 

2. Gluons can carry a very large part of the nucleon spin. If we assume that 

1 

J Ag(x, Q 2 = W)dx « 2 (2.42) 
o 

then [13]: 

l 

J As(x, Q 2 )dx = -0.06 ± 0.06 (2.43) 
o 

l 

J AS(x, Q 2 )dx = 0.42 ± 0.17 (2.44) 
o 

3. The polarised sea quarks distribution function is large and negative, so 
it cancels most of the valence quark contribution to the nucleon spin. 

4. The orbital momenta of the nucleon constituents L Z (Q 2 ) is large and it 
is a lacking part of the nucleon spin, according to the total nucleon spin law 
conservation: 

i(AS) + (A 5 )+L z (Q 2 ) = i (2.45) 

The problem of the "spin crisis" is still open. Its solution requires most of 
all the knowledge about polarised distribution functions of gluons. Theo- 
retical analyses QCD as well as future polarised experiments at HERA [8], 
COMPASS [15] and RHIC [21] should better determine the gluon spin func- 
tion behaviour, particularly in the small Bjorken x region, what is a crucial 
point to overcome the " spin crisis" . 

3. Spin structure functions in the small Bjorken x region 

Determination of the nucleon spin structure functions in the small Bjor- 
ken x region is very important from both theoretical and experimental point 
of view. Because of technical limit, present experiments do not give any in- 
formation about small x region (x ~ 10~ 4 , 10~ 5 ) and therefore there are 
still the uncertainties in the determination of parton distribution functions 
(in particular the gluons) in this region. Theoretical analyses, based on the 
perturbative QCD, allow to calculate the nucleon structure function within 
some approximations (Q 2 LO, Q 2 NLO, lnl/x etc.). Choice of some particu- 
lar approximation depends of course on the region of its application and the 
basic criterion is agreement of the theoretical predictions with experimental 
data. Thus the small x behaviour of the nucleon spin structure functions 
implied by QCD can be via sum rules (BSR, EJSR) tested experimentally. 
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Moreover, the aim of the QCD analysis is to yield an adequate, compact 
description of the nucleon structure functions in the whole range of x: for 
the values of x, which are available in experiment and which are not as well. 
The small x region is also a challenge for QCD analysis, because theoretical 
predictions of the structure function g\(x, Q 2 ) at low x are relevant for the 
future polarised HERA measurements [8]. 

Structure functions for small x and fixed Q 2 are connected with the 
virtual Compton scattering total cross-section at very high energy W 2 — > oo: 

W 2 = Q 2 {- - 1) (3.1) 

x 

where W is the total CM energy. The small value of x (x — ► 0) corresponds 
by definition to the Regge limit and therefore the small x behaviour of 
structure functions can be described using the Regge pole exchange model 
[2]. The Regge theory predicts, that spin dependent structure functions 
gP,n,d - n sma n x region behave as 

g p,n,d ^ x - a (3.2) 

where a denotes the axial vector meson trajectory and lies in limits: 

-0.5 < a < (3.3) 

The experimental data from HERA confirm such a Regge behaviour of struc- 
ture functions (3.2) but only in the low Q 2 region Q 2 < A 2 (A 2 « 200 MeV) 
i.e. in the region, where the perturbative methods are not applicable. At 
larger Q 2 , because of parton interaction, the structure functions undergo 
the Q 2 evolution and their behaviour, implied by perturbative QCD is more 
singular than that, coming from the Regge picture. This fact is also in agree- 
ment with experiments of unpolarised as well as polarised DIS. It is well 
known at present, that for x — > the Regge behaviour x~ a (—0.5 < a < 0) 
is less singular than the perturbative QCD predictions for all of parton dis- 
tributions except unpolarised, nonsinglet (valence) quarks qNS- 
There are two basic approaches of perturbative analysis of unpolarised struc- 
ture functions in the small x region: a resummation of the GLAP Q 2 loga- 
rithms [3], [10], [16]: 

£MQ 2 )ninQ 2 r (3.4) 

n,m 

or a resummation of 1/x logarithms BFKL [4]: 

]TMQ 2 )ninir (3-5) 

n,m ^ 
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In the latest experimental data analyses of DIS the standard approach is 
based on the QCD evolution equations within next-to leading (NLO) ap- 
proximation [3], [10]. This approach is appropriate for polarised as well as 
unpolarised nucleon structure functions. It has to be emphasised, that the 
changes, which appear by the transition from the LO approach to the NLO 
one are different in the case of polarised and unpolarised structure func- 
tions. Splitting functions P^\x) and AP^\x), governing the evolution 
of quark and gluon distribution functions in the LO approximation have a 
form (A.7)-(A.10) in the spin independent case and a form (A.11)-(A.14) in 
the spin dependent case respectively. In the small x limit x — ► they take 
a LO form: 

P<£>(*)~!; n ( §(*)-4 P G>)~h P GG(*)~P (3-6) 
APW(x)~|; AP$0r)~-i; APg(x) ~ |; AP® (x) ~ 9 (3.7) 

In the NLO approximation, functions P^p(x) and AP^\x) have the fol- 
lowing small x behaviour [16], [17]: 

^(^)~~; (3-8) 

J/ J/ tiy tiy 

AP$(x), AP$(x), APg(x), ApW(x)~a + 61nx + cln 2 x (3.9) 

where a,6,c are constants. Prom the comparison of (3.6)-(3.9) one can read 
an interesting information about the small x behaviour of polarised and 
unpolarised nucleon structure functions. First: comparing (3.6) with (3.8) 
i.e. spin independent Pij(x) function in LO and NLO approximations, one 

can notice, that the singular terms 1/x appear in Pqq\x) and P^q(x) while 

they do not in Pqq (x) and P^q (x) respectively. This means that the evolu- 
tion of unpolarised, fermion distribution functions at small x is completely 
dominated by the NLO terms in the perturbative calculation. Secondly: 
the singular l/x terms in Pq^{x) and Pqq(x) are opposite in sign with re- 
spect to the corresponding LO terms: Pq^(x) and Pqq{x). Thus, the rapid 
growth of the unpolarised gluon distribution function at small x, well known 
from LO analyses, will be dampened by the NLO contribution. Thirdly: all 
spin dependent splitting functions AP-j\x) are nonsingular (without l/x 

terms). (An exact form of AP-j\x) is given e.g. in [16].) It means, in 
general, that in the small x region spin dependent parton distribution func- 
tions are by a factor of x less singular than the corresponding unpolarised 
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distributions. Knowledge of the nucleon structure function at small x is 
extremely important: gi(x,Q 2 ) results for this region enter into integrals 
for moments T\ ,n,d (2.10) and hence into the Bjorken and Ellis-Jaffe sum 
rules (2.20)-(2.23). From the other side, lack of the very small x (x < 10 -4 ) 
experimental data causes that the all knowledge of structure functions in 
limit small x (x — > 0) comes almost entirely from theoretical QCD analyses. 
Therefore perturbative (and nonperturbative as well) QCD predictions in 
the small x region are at present of great importance. 

Evolution equations GLAP for unpolarised gluon distribution functions 
g(x,Q 2 ) and unpolarised singlet quark distributions T>(x,Q 2 ) 

N f 

E(ar, Q 2 ) = 5>i(z, Q 2 ) + %(*, Q 2 )] (3.10) 
i=i 

lead to the following behaviour [18]: 

n 

xg(x, Q 2 ) ~ a' l / 2 e 2 ^-^(l + ]T eV+Vj (3.11) 

i=i 

n 

xT(x,Q 2 ) ~ p- x a- x l 2 e 2 ^-^(\ + Y J t l fp i+l <) (3.12) 

i=i 

where £, £, P are given as: 

n = in the summation (3.11) and (3.12) corresponds to LO while n = 1 to 
NLO approximations. Remaining parton distributions i.e. p = qNSi ^QNSj 
AS, Ag behaves in the small x region like: 

n 

p{x, Q 2 ) ~ ff- 1 /2 e 27/^/C( 1 + ^ 4p 2i+1 a l j (3.14) 

i=l 

One can read from (3.11), (3.12) and (3.14), that spin dependent parton 
distribution functions xAqjys, xAS, xAg and the unpolarised nonsinglet 
quark distribution xq^s are by a factor of x less singular than the unpo- 
larised distribution xg and xS. Moreover, in the case of p function (p = qNS, 
Aqws, AT,, Ag) (3.14) higher order corrections are more important than 
in the xg (3.11) and (3.12) case because of appearance p 2l+l terms in 
(3.14) instead of p t+1 in (3.11) and (3.12). It must be also emphasised, 
that the small x behaviour of parton distributions strongly depends on the 
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input parametrisation p(x, Qq) (2.32). When it is nonsingular, then the 
singular small x behaviour of parton distributions is fully generated by Q 2 
evolution and has a form (3.11)-(3.14). Whereas in a case of singular input 
parametrisation e.g. [7] 

Ag(x,Q 2 )^x a (l-xf(l + 7 x 5 ) (3.15) 

where a = —0.5, (3 = 4, 7 = 3, <5 = 1, singular x a small x behaviour of Ag 
distribution will survive Q 2 QCD evolution and will be leading towards the 
singular terms, generated perturbatively. Choosing the input parametrisa- 
tion (2.32), experimental data must be taken into account. As a start scale 
for the GLAP evolution equations Qq = 1 GeV 2 is assumed at present. This 
value is a limit where from the one side (Q 2 < Qq) experimental measure- 
ments confirm the Regge theory of DIS, and from the other side Q 2 > Qq it 
is a starting point for methods of perturbative QCD (Q 2 > 1 GeV 2 ). There 
are two approaches dependent on the choice of the input parametrisation. 
Either we assume Regge behaviour (3.2) of parton distributions at small x 
and then their shape at small x and larger Q 2 Q 2 > Qq is fully implied 
by NLO QCD evolution (3.11)-(3.14) or we take singular input parametri- 
sations (but no more singular than it results from experimental data i.e. 
AqNs( x ,Q 2 ) < x -0 ' 5 ), which will survive LO and NLO GLAP evolution. 
The Regge theory predicts the following behaviour of parton distributions 
at small x and Q 2 < 1 GeV 2 : 

x£ ~ const (Pomeron) 
Qns ~ x -0 - 5 (Reggeon A 2 : p - u)\ 
AE,Aq NS ~ x° -hx - 5 (Reggeon A x ) (3.16) 

Experimental analyses [5], [19] are in a good agreement with (3.16) for 
Q 2 < 1 GeV 2 and also for larger Q 2 , after taking into account the perturba- 
tive effects. Investigating small x region within perturbative methods, one 
should include all of those terms in C(x, Q 2 ) and Pij(x, a s ), which remark- 
ably influence the shape of the nucleon structure functions and what will 
be soon verified experimentally. It has been lately noticed [20], [24], [25] that 
the spin dependent structure function g\ in the small x region is dominated 
by ln 2 (l/x) terms. These contributions correspond to the ladder diagrams 
with quark and gluon exchanges along the ladder - cf Fig. 2. The contri- 
bution of non-ladder diagrams to the nonsinglet spin dependent structure 
function is negligible. Thus the behaviour of the spin dependent nucleon 
structure functions at small x is expected to be governed by leading double 
logarithmic terms of type a"ln 2n (x). These terms must be resummed in 
the coefficients and splitting functions Pij(x,a 2 ). The resummation of the 
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Fig. 2. A ladder diagram generating double logarithmic ln 2 (l/x) terms in the non- 
singled spin structure function g\. 



double logarithmic terms In 2 x in the limit of a very small x (x — ► 0) is given 
by the following equations [20]: 

1 k 2 /z 

fNs(x,k 2 ) = ff s (x,k 2 ) + ^-j^ J ^-AP£\z)f NS (- z ,k' 2 ) (3.17) 



X 1,2 
K 



1 k 2 /z 

a, f dz f dk' 2 



X 1,2 

Ap(y(z) /s (^' 2 )+AP ? 5( 2 )//,fc' 2 ) 



(3.18) 



1 k 2 /z 

a, f dz f dk' 2 



fg (x,k 2 ) = fP(x,k 2 ) + ^-Jf J 

X 1,2 

ap(J( z )/ s (^' 2 ) + apS(z)//,^; 



(3.19) 



17 



where AP^\z) are LO approximation at z —> and have a form (3.7). 
Unintegrated distributions / in equations (3.17)-(3.19) are related to the 
corresponding polarised distributions Ap(x,Q 2 ) via 

W 2 2 

Ap(x,Q 2 ) = ApW(x) + J ^f(x' = x(l + k 2 /Q 2 ),k 2 ) (3.20) 

k 2 



c Ah 2 

Ap(°)(x) = y ^/(x,fc 2 ) (3.21) 
o 



A; 2 is the transverse momentum squared of the parton. The inhomogeneous 
term /(°)(a;,fc 2 ) in the nonsinglet case (3.17) has a form 

r (0) / l.2x _ «*(^) 4 + ^ 2 ) Ap (0) (x/z) - 2zAp(°) (x) 



f NS (x, k ) 



2f 3J z 1 — 2 



+ ^(M [2 + |i n (i_ x )]Ap(°)(x) (3.22) 
2ir 3 



where the nonperturbative part of the distribution Ap^ (x) is parametrised 
on the basis of the experimental data at the small values of Q 2 (Q 2 < 
1 GeV 2 ). This nonperturbative parametrisation is given by 

A P P s (x) = N(l - xf (3.23) 

Parameter 77 in the valence quark case (Aq^s) is equal to 3, gA = 1.257 
is the axial vector coupling and N is the normalisation constant, in the 
nonsinglet case determined from the Bjorken sum rule (2.11), (2.12), which 
can be written as: 

l l 
J g? s (x,Q 2 )dx = J{g p l -g r l)(x,Ql)dx = ^g A (3.24) 
oo 

The source of the double logarithmic terms In 2 x in gi(x, Q 2 ) is the double 
integration in the formula for function f(x,k 2 ): 

X 1.2 
K Q 
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where the upper limit in the integral over the transverse momentum k 
is z-dependent (= k 2 jz). Thus, double logarithmic terms come from the 
integration over the longitudinal momentum fraction z together with the 
integration over k' 2 with z-dependent upper limit: 

f(x, k 2 ) ~ ln 2 (l/x) = In 2 x (3.26) 

Equations (3.17)-(3.19) in the case of the fixed coupling constant a s can be 
solved analytically [20]. These equations generate singular small x behaviour 
of the polarised parton distributions and hence of the spin dependent struc- 
ture function g\ i.e.: 

g? s (x,Q 2 )~x- x » s ; gf(x,Q 2 ) ~ x~ Xs ; Ag(x, Q 2 ) ~ x~ Xs (3.27) 

where g^ s = g{ — g™ , gf = g\ + g™ , Ag is the polarised gluon distribution 
and exponents Aj have forms: 

Xns = 2 ] I^AP^\x); A s = 2^7 (3-28) 

7 + = I[APW(x) + ApW(x) 

+ 0AP4 o) (x) - APg^P + 4AP5(x)APg(x)] (3.29) 

As it has already been mentioned above, the singular small x behaviour of 
the polarised structure function (3.27) become leading only in the case of the 
nonsingular input parametrisation (2.32) e.g. for the simple parametrisation 
(3.23). Because only unpolarised, nonsinglet parton distributions qjys have 
Regge small x behaviour a; -0 - 5 more singular than that, implied by QCD 
evolution, the shape of all spin dependent distributions is mostly governed 
by QCD evolution. Thus the leading small x behaviour of polarised nucleon 
structure functions is (3.27) and moreover singlet structure functions dom- 
inate over the nonsinglet ones (As > Xns)- For example, at Q 2 = 10 GeV 2 , 
Nf = 3 and scale parameter A = 232Mey, the corresponding values of Aj 
are: = 1.4, Xns = 0.48 (a s = 0.27, 7 + = 11.7). It must be pointed out, 
that the equations (3.17)-(3.19) can be applied only for very small x. They 
are based on the approximation, that for very small x splitting functions 
AP^j\x — ► 0) are given by (3.7). At large and moderately small values of x 
this approach is no more adequate. In the region of "not too small x" , which 
can be explored experimentally, in theoretical QCD analyses one should use 
the GLAP (LO or NLO) equations with complete APij(z) functions. Com- 
bining the standard LO GLAP approach with the double In 2 x resummation, 
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it is possible on the one hand to guarantee an agreement of QCD predictions 
with experimental data in the large and moderately small x and on the other 
hand to generate the singular small x shape of polarised structure functions, 
governed by In 2 a; terms. With this aim the equations (3.17)-(3.19) should 
be extended to include complete splitting functions APij(z) and not only 
their approximations (3.7) for z — > 0. In this way one can obtain system of 
equations, containing both LO GLAP evolution and the double logarithmic 
In 2 x effects at small x. Analyses of such unified GLAP LO + In 2 x approach 
are presented in [20] . On the basis of this interesting method we give in the 
next chapter the predictions for the gi S function in the case of nonsingular 
as well as singular input parametrisation Apns(x, Qq). 



4. Predictions for the nonsinglet spin structure function g± 

The small x behaviour of both nonsinglet and singlet spin dependent 
structure functions g^ s (x,Q 2 ) and gf(x,Q 2 ) is governed by the double 
logarithmic terms a™ln 2n (x) [20], [24], [25]. But in contrast to the singlet 
polarised function, for the nonsinglet one the contribution of nonladder di- 
agrams is negligible. Thus we should consider only ladder diagrams with 
quark (antiquark) exchange, Fig. 2. Hence the nonsinglet part of the po- 
larised structure function gi has a form: 

g? s (x,Q 2 )=g p 1 (x,Q 2 )-g?(x,Q 2 ) (4.1) 

where g\ and g™ are spin dependent structure functions of proton and neu- 
tron respectively. According to (2.3), (2.4), (2.29)-(2.31) one can obtain for 
the colour number N c = 3: 

gPl = \ Au + Ts Ad + Is A ^ + l A * (42) 
T8 Au + l Ad + T8 Au + l A " (43) 

i ( Au val - Ad val ) = ^ (An - Ad) (4.4) 
(4.4) results from the assumption of SU(3) flavour 

Au = Ad (4.5) 

and hence all of gluon and sea quark contributions from the proton and the 
neutron structure function cancel mutually. This feature that the small x 
behaviour of the spin dependent nonsinglet structure function is governed by 
the double logarithmic terms a n s ln 2n (x) is very important from the point of 



and hence 



9 n i 



9? S 



The simple form of g± 
symmetry: 
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view of small x QCD analysis. This is different from the case of unpolarised 
nonsinglet structure functions F^ 5 , where the small x behaviour of F2, 
generated by the a™ln 2n (x) terms, is dominated by the nonperturbative 
contribution of A2 Regge pole. For gi S the relevant Ai Regge pole has low 
intercept a7vs(0) < and for small x in the Regge limit one has: 



Thus the Regge behaviour of the spin dependent structure functions is un- 
stable against the resummation of the In 2 x terms, which generate more 
singular x shape than that (4.6) with ajvsCO) < 0. Therefore the mea- 
surement of the nonsinglet spin dependent structure function can be a very 
important test of the QCD perturbative analyses in the small x region. In 
our numerical analysis we follow [20] and [25]. Solving the unified equa- 
tion incorporating GLAP Q 2 evolution and the In 2 x resummation we get 
the results for the nonsinglet polarised structure function (x, Q 2 ) in the 
perturbative region Q 2 > Qq for different values of x G (0; 1). This equation 
taking into account both GLAP evolution and In 2 x effects for g^ s function 
has a form [20], [25]: 



9? S (x,Q 2 ) 



x 



ajvs(O) 



(4.6) 




k 2 /z 




+ ( 1 - + ^ln(l-x))f(x,k' 2 )} 



(4.7) 



where 




X 



+ (- + £ ln(l - x))sH*)] 



(4.8) 



g 2 (i/x-i) 





(4.9) 



and 




.2 



(4.10) 







21 



Comparing (4.7)-(4.10) with (3.17), (3.20)-(3.22) it is clear that AP qq (x) in 
(4.7) has a full GLAP form instead of its approximation for x — > in (3.17) 
and g^ s plays simply the role of Ap from (3.20). We solve eq.(4.7) using 

different parametrisations of g^ S ^°\x): the simple one, implied by Regge 
behaviour of g^ s in nonperturbative region 

9? S{0 \x) = 9? S (x,Ql)=N(l-xf (4.11) 

and two dynamical inputs: GRSV (Gliick, Reya, Stratmann, Vogelsang) [22] 
and GS (Gehrmann, Stirling) [23]. The nonsinglet spin dependent structure 
function must satisfy the Bjorken sum rule (2.11), (2.12) independently of 
the value of Q 2 . This means that for any Q 2 , the first moment of gi S must 
be equal to l/6gA similarly to the case of the low scale Qq (3.24): 

i l 
(9? S (*, Q 2 )) = J 9? S (x, Q 2 )dx = /(<£- 9?){x, Q 2 )dx = l -g A = 0.2095 



(4.12) 

This condition implies the proper normalisation constants N in all of input 
parametrisations. Thus the constant TV in (4.11), found from the Bjorken 
sum rule is equal to 2/3gA = 0.838 (we set gA = 1.257) and the Regge 
nonsingular input (4.11) takes a form: 

REGGE : g? s (x, Q 2 ) = \g A {l - xf = 0.838(1 - xf (4.13) 

The Regge behaviour of structure functions at small x, as it was men- 
tioned above, has been confirmed by HERA experiments in the low Q 2 
region (Q 2 < 1 GeV 2 ). Therefore the choice of the Regge input allows to 
unite the nonperturbative origin with QCD perturbative analysis starting 
at Qq ~ 1 GeV 2 . In this way, assuming the Regge (flat, nonsingular) be- 
haviour of structure functions at low Q 2 scale i.e. Qq = 1 GeV 2 , we expect 
that the singular small x behaviour of polarised structure functions is com- 
pletely generated by QCD evolution, involving NLO or even (as in our case) 
GLAP+ln 2 x approach. This analysis, based on the Regge input (4.13), is 
however one of two main possible scenarios, describing the small x behaviour 
of spin structure functions. The second is to allow steeper (more singular) 
inputs of structure functions at Qq, what intensifies more the growth of 
structure functions as x — > implied by QCD. The only constraint on these 
two scenarios is consistency of their predictions with experimental data. In 
our analysis of the gi S structure function we consider dynamical inputs pro- 
posed by GRSV [22] and GS [23]. These inputs result from a global analysis 
of all available recently deep inelastic polarised structure function data [14]. 
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Our calculations incorporating both GLAP evolution and resummation of 
the In 2 x terms are based on the LO fitted inputs. In such a way the spin 
dependent nonsinglet structure function g^ s (4.4) has an input form: 

GRSV : g? s {x, Q 2 = 1 GeV 2 ) = 0.327x^ a267 (l - 0.583x a175 + 1.723z 
+ 3.436x 3 / 2 )(l - x) 3 - 486 + 0.027x- a624 (l + 1.195x a529 
+ 6.164x + 2.726x 3/2 )(l-x) 4 - 215 (4.14) 

GS : g? S (x, Q 2 = 4 GeV 2 ) = 0.29x" a422 (l + 9.38x - 4.26^) 

x (l-x) 3 - 73 + 0.196x~ a334 (l + 10.46x- 5.10^(1 -x) 4 ' 73 (4.15) 

for details see Appendix B. All numerical calculations have been performed 
in C code on PC computer under LINUX system. Our numerical re- 
sults for g? s based on Regge (4.13), GRSV (4.14) and GS (4.15) input 
parametrisations are presented in Figs. 3-7. In Fig. 3 we plot different input 
parametrisations (x, Qq). Figs. 4,5 show the nonsinglet function gi' s af- 
ter evolution to Q 2 = 10 GeV 2 for these different parametrisations (Regge, 
GRSV, GS) and Figs. 6, 7 present the function 6xg^ s = x(Au va i — Ad va i) at 
Q 2 = 10 GeV 2 also for different inputs gi S (x, Qq). In all of Figs. 4-7 pure 
GLAP evolution is compared with double logarithmic In 2 x effects at small 
x. Additionally, in Figs. 5-7 we compare our numerical results with recent 
SMC (1997) data [14]. Contributions 6(g? s ) (4.12) and 6AI(x a ,x b ,Q 2 ) 

x b 

AI(x a , x b , Q 2 ) = J g? s (x, Q 2 )dx (4.16) 

Xa 

to the Bjorken sum rule at Q 2 = 10 GeV 2 together with experimental SMC 
values are presented in Tab. 2: *) means the extrapolation of experimental 
data to low x and **) is the integral over the measured range of x. From 
Figs. 4-7 one can read that the double logarithmic In 2 x effects are very sig- 
nificant for x < 10~ 2 . Besides, as it has been expected, the growth of the 
nonsinglet proton spin structure function g? s in the very small x region is 
much steeper for dynamical parametrisations (GRSV or GS) than for the 
Regge one. The comparison of our theoretical model with experimental data 
in Tab. 2 and Figs. 5-7 yields the conclusion that all of the theoretical predic- 
tions for different parametrisations (Regge, GRSV, GS) and incorporating 
pure LO GLAP QCD evolution as well as LO GLAP evolution with In 2 x 
effects are in a good agreement with experimental data within statistical 
errors. Unfortunately, the most interesting x region is still nonavailable for 
experiment. So the problem, which QCD approach is the most adequate for 
the description of small x physics in the polarised deep-inelastic scattering 
of particles remains unsolved. 
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TABLE II 



PARAMETRISATION 


6AI 


6A7 


6A7 




(0,1, Q 2 ) 


(0, 0.003, Q 2 ) 


(0.003,0.7, Q 2 ) 


INPUT 


1.257 


0.0150 


1.232 


REGGE LO GLAP 


1.255 


0.0342 


1.219 


LO GLAP+ln 2 x 


1.249 


0.0493 


1.198 


INPUT 


1.257 


0.0786 


1.194 


GRSV LO GLAP 


1.249 


0.107 


1.171 


LO GLAP+ln 2 x 


1.242 


0.119 


1.153 


INPUT 


1.257 


0.123 


1.160 


GS LO GLAP 


1.253 


0.134 


1.151 


LO GLAP+ln 2 x 


1.247 


0.142 


1.139 


EXPERIMENT 


1.29+0.24 


*) 0.09+0.09 


**) 1.20+0.24 



Table 2. Theoretical contributions 6A7(x a , Xb, Q 2 ) and their experimental SMC 
values 
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Fig. 3. Input parametrisations of the nonsinglcd spin structure function of the 
proton g? s (x,Ql): REGEE (4.13) - dotted line; GRSV (4.14) - solid line; GS 
(4.15) - dashed line. 
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Fig. 4. 5 f 5 at Q 2 = 10 GeV 2 based on inputs: REGGE - dotted, GRSV - solid, 
GS - dashed. For each pair of lines the pure LO GLAP prediction lies below the 
LO GLAP+ln 2 x one. 
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Fig. 5. g^ s at Q 2 = 10 GeV 2 ; similarly as in Fig. 4 but for the measurable experi- 
mentally region of x. Squares show the recent SMC data 1997 [14]. 
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Fig. 6. Function 6xgf s at Q 2 = 10 GcV 2 . Predictions based on the REGGE input 
- dotted and GRSV - solid. LO GLAP above LO GLAP+ln 2 x at x = 0.2. SMC 
1997 data with statistical errors are shown. 
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Fig. 7. LO GLAP+ln 2 x predictions for function 6xg^ s at Q 2 = 10 GeV 2 , based 
on input parametrisations: REGGE - dotted, GRSV - solid, GS - dashed. Plots 
are compared with SMC 1997 data. 
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5. Summary and conclusions 



In this paper the main theoretical and experimental problems in nucleon 
spin structure physics have been briefly reviewed. The results of current 
experiments are deviation of the Ellis-Jaffe sum rule and validity of the 
Bjorken sum rule. This causes that the question "how is the spin of the 
nucleon made out of partons?" is still open. The great puzzle are exper- 
imental results which violating the Ellis-Jaffe sum rule imply that only a 
very small part of the spin of the proton is carried by quarks. So where is 
the nucleon spin? Maybe gluons take a large fraction of the nucleon spin? 
Or maybe the spin of the proton is "hidden" in orbital angular momentum 
of quarks and gluons? Maybe at last the solution of the spin problem lies 
in the small x physics and the lacking spin of the nucleon is hidden in the 
unmeasured very small x region. The answer the above questions will be 
possible thanks to the progress in theoretical and experimental research in 
the small x physics. Perturbative QCD analysis, based on GLAP evolution 
equations is in a good agreement with experimental data. This agreement 
concern unpolarised and polarised structure functions of the nucleon F\, 
F2, g\ within NLO approximation in the large and the moderately small 
Bjorken x region. Unfortunately, practically lack of the experimental mea- 
surements in the very small x region (x < 10~ 3 ) makes the satisfactory 
verification of the theoretical QCD predictions in this region impossible. 
Knowledge of the behaviour of the nucleon spin structure functions when 
x — > is crucial in the determination of Bjorken and Ellis-Jaffe sum rules 
i.e. in overcoming the "spin crisis". Understanding of the small x physics 
in the polarised DIS processes requires to take into account all of these 
perturbative QCD effects which become significant in the small x region 
and which could be verified by future experiments. Present QCD analy- 
ses, based on the GLAP LO or NLO Q 2 evolution seem to be incomplete 
when x — ► 0. The growth of the unpolarised as well polarised structure 
functions of the nucleon in the small x region is governed by leading dou- 
ble logarithmic terms of the form a™ ln 2n (x) , generated by ladder diagrams 
with quark and gluon exchange. This singular behaviour of the structure 
functions at low x, implied by In 2 x terms, is however better visible in the 
polarised case. For unpolarised, nonsinglet structure functions of the nu- 
cleon the QCD evolution behaviour at small x is screened by the leading 
Regge contribution. Therefore the spin dependent structure functions of the 
nucleon are a sensitive test of the perturbative QCD analyses in the low x 
region. Our numerical analyses incorporating the LO GLAP evolution and 
the In 2 x effects at small x show, that the growth of the nonsinglet polarised 
structure function of the nucleon g^ s , implied by In 2 x terms, is significant 
for x < 10 -2 . Our predictions for g^ s are in a good agreement with the 
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recent SMC data for small x region (x ~ 10~ 3 ). The contribution from the 
low x region (x < 0.003) to the Bjorken sum rule is found to be around 4% 
(for Regge input g^ s (x, Qq)) and 10% (for dynamical inputs) of the value 
of the sum. Theoretical predictions for gi' n,d , taking into account the In 2 x 
resummation effects will be in the future verified experimentally. There are 
a few hopeful experimental projects of the investigation of the nucleon's spin 
structure. One of these is the HERMES experiment (start in 1995) located 
in HERA at DESY with a fixed polarised H,D or 3 He target and a longitu- 
dinally polarised positron beam of 27.5 GeV [26]. The accessible kinematic 
range is 0.004 < x < 1 and 0.2 < Q 2 < 20 GeV 2 . The HERMES experiment 
allows a direct measurement of the polarised quark distributions for individ- 
ual flavours also g\ ,n,d {x, Q 2 ) and even g2(x, Q 2 ). The question of the gluon 
polarisation is also addressed experimentally. The polarised gluon distri- 
bution Ag(x,Q 2 ) may play a crucial role in understanding of the nucleon 
spin structure. The measurement of Ag(x, Q 2 ) in the charm production via 
photon-gluon fusion process j*g — > cc will be possible at COMPASS exper- 
iment at CERN [15]. In this project the polarised muons will be scattered 
on polarised proton and deuteron targets. The energy of the muon beam 
will be of 100 GeV and 200 GeV and the Bjorken x region x > 0.02. The 
COMPASS measurements are expected to start in 2000. A very important 
program which will test many elements of QCD in the perturbative as well 
as in the nonperturbative region is RHIC spin project at Brookhaven [21]. 
This program with polarised proton-proton collider will start in 2000 and 
will allow for a measurement of the polarised gluon density via heavy quark 
production (gg — > QQ) or via direct photon production {gq — > ^yq). Finally, 
a very promising experimental project in high energy spin physics is HERA 
[8]. The polarisation of the proton and electron beams at y/s = 300 GeV 
will enable to measure the structure function g\(x, Q 2 ) and spin dependent 
quark distributions Aqf(x,Q 2 ) at very low x (x ~ 10~ 5 ). From polarised 
di-jet production it will be possible to determinate the polarised gluon dis- 
tribution Ag(x, Q 2 ) for the region 0.002 < x < 0.2. Additionally in HERA, 
a program of polarised proton-proton collisions is proposed. This high en- 
ergy proton-proton scattering will allow via J ftp production for the direct 
determination of the gluon function Ag(x,Q 2 ). The new HERA with po- 
larised experiments and the largely extended kinematical region of x and 
Q 2 will contribute a lot to our understanding of high energy spin physics. 
The problem of the spin structure of the nucleon is nowadays one of the 
most important challenges for theory and experiment as well. 
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Appendix A 

GLAP evolution of polarised quark and gluon distribution functions in the 

nucleon 

In perturbative QCD distribution functions of partons q(x, Q 2 ), g(x, Q 2 ), 
Aq{x, Q 2 ), Ag(x, Q 2 ) evolve with Q 2 . This evolution is described by GLAP 
[3], [9] equations and it is assumed that for polarised functions Ap(x, Q 2 ) 
the evolution equations have the same form like for unpolarised functions 

p{x,Q 2 ): 

j t [Aq NS (x,t)} = ^-AP qq ®Aq NS (x,t) (A.l) 
±[Aq s (x,t)] = ^-[AP qq ®Aq s (x,t)+2N f AP qG ®Ag(x,t)} (A.2) 

at Z7T 

^[Ag(x,t)] = ^[AP G5 Aqs(x,t) + 2N f AP GG ® Ag(x,t)} (A.3) 



where t = ln(Q 2 /A 2 ) and 

P(x)®q(x,t)= [ ^P(z)q(-,t) (A.4) 

J z z 



1 

dz 



Functions C(x,a s ) and Pij(x,a s ) are calculated in leading LO approxima- 
tion or with the next correction to LO i.e. in NLO approximation with 
respect to coupling constant a s : 

LO: P t3 (x,a s ) = a s P t f(x) (A.5) 

NLO : Pij{x, a s ) = a s P^\x) + a 2 a Pjp(x) (A.6) 

Functions Pij(x) for the unpolarised case are different from those APij(x) 
for the polarised one [10]. In LO they have a form: 

p « )(x) = l(T^: + 26(1 ~ x) (A - 7) 

P${x) = \{x 2 + {l-x) 2 ) (A.8) 
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P ^ W = 3[ (T^ + V^^^ + r 50 "' ;! ,A ;ni 

and 



"W=l<rk (A - n) 

AP 9 ( g(x) = l(2x-l) (A.12) 



AP^(x) = -(2-x) (A.13) 

AP ^ {X) = 3[ (T^ + (3 " 3x + x 2 + x 3 ) - Ls(l - x)} (A.14) 
where (1 — x) + is denned as: 

1 /(*)<;* _ l/M-Zd)^ 



4 



o (i-x)+ y 



Appendix B 



Dynamical input parametrisations of the nonsinglet polarised structure 

function g^ s 

In our calculations we adopt GRSV (Gliick, Reya, Stratmann, Vogel- 
sang) [22] and GS (Gehrmann, Stirling) [23] parametrisations of valence 
quarks Aq va i. We assume SU(3) flavour symmetric scenario, where 

Au = Ad (B.l) 

This assumption leads to formula (4.4): 

9i S = g{-g r l = \{Au- Ad) = ^(Au val - Ad val ) (B.2) 

Input parametrisation of Au va i and Ad va i have a general form: 

GRSV: Aq va i = Nx a2 x ai - 1 {l + Ax b + Bx + Cx 3 / 2 ){l-x) D (B.3) 
GS: Aq val =N'x a '- 1 {l + 1 x + p^){l-x) D ' (B.4) 
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where N, N' are normalisation factors, implied by Bjorken and Ellis-Jaffe 
sum rules (2.11)-(2.14). These sum rules for input scale Qq can be read as 

1 

a 3 = J {Au val - Ad val )dx = 1.257 (B.5) 
o 

l 

a 8 = J (Au va i + Ad va i)dx = 0.579 (B.6) 
o 

(B.5) and (B.6) give immediately 

l 

J Au val dx = 0.918 (B.7) 
o 

l 

J Ad vai dx = -0.339 (B.8) 
o 

what allows to find N, N' factors. The full set of input parameters for 

GRSV and GS distributions is as follows: 

GRSV: 

Ql = 1 GeV 2 , Aqcd = 232 MeV 

for Au val : N = 1.964, a x = 0.573, a 2 = 0.16, b = 0.175, 

A = -0.583, B = 1.723, C = 3.436, D = 3.486 
for Ad val : N = -0.162, ai = 0.376, a 2 = 0, b = 0.529, 
A = 1.195, B = 6.164, C = 2.726, D = 4.215 

GS: 

Q 2 = 4 GeV 2 , Aqcd = 200 MeV 

for Au val : N' = 1.741, a' = 0.578, 7 = 9.38, p = -4.26, D' = 3.73 
for Ad vat : N' = -1.176, a' = 0.666, 7 = 10.46, p = -5.10, D' = 4.73 

In both GRSV and GS inputs we employ the LO fits. Thus the input 
parametrisations have final forms: 

GRSV: 

Au va i = 1.964x-°- 267 (l-0.583x ai75 + 1.723x + 3.436x 3/2 )(l-x) 3 - 486 (B.9) 
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(B.10) 



Ad va i = -0.162x~ a624 (l + 1.195x a529 + 6.164x + 2.726x 3/2 )(l 
GS: 

Au va i = 1.741x" - 422 (l + 9.38x - 4.26y/x)(l - xf' 73 (B.ll) 

Ad val = -1.176x"°- 334 (l + 10.46x - 5.10>/x)(l - x) 4 ' 73 (B.12) 
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